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We present and discuss a list of some interesting points that are currently open in nonextensive 
statistical mechanics. Their analytical, numerical, experimental or observational advancement 
would naturally be very welcome. 
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q ! 1. Introduction 

. Statistical mechanics is one of the pillars of contemporary physics. Indeed, it provides in principle the 
r* \ link between the microscopic and the macroscopic description and understanding of the world through a 
. — , ■ rich variety of mesoscopic instances. More precisely, starting from mechanics (classical, quantum, relativis- 
^ . tic) and electromagnetism, it incorporates appropriate concepts of the theory of probabilities, and finally 
leads, for large enough systems, to thermodynamics. It does so through various intermediate descriptions 
such as those involving master, Langevin, and Fokker-Planck equations. See Fig. [1] for a schematic set of 
connections. 

The statistical mechanical approach of physical phenomena was first introduced by Maxwell, Boltz- 
mann and Gibbs. A crucial step was done by Boltzmann who proposed a specific functional form for the 
thermodynamic entropy (that had been introduced a few years earlier by Clausius) in terms of the proba- 
bilities of the microscopic configurations. Its expression for say a system describable through a discrete set 
of probabilities {pi} of W possibilities is given by 

w w 

Sbg = -k^pilnpi (^ P j = l), (1) 

i=l i=l 

where BG stands for Boltzmann- Gibbs, and k is a positive constant referred to as the Boltzmann constant 
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Fig. 1. Relevant connections between microscopic (bottom; red boxes) and macroscopic (top; blue box) descriptions through 
mesoscopic (middle; orange, green and black boxes) ones. 



(which, together with the velocity of light c, the Planck constant h and the gravitational constant G, 
constitutes a minimal set of universal physical constants); for convenience it is sometimes taken k = 1. The 
important particular instance where probabilities are all equal (i.e., pi = 1/W, Vi) leads to 



Sbg = khiW , 



(2) 



carved on Boltzmann's grave stone in the Central Cemetery in Vienna. It is straightforward to verify that 
in all cases Sbg > 0. 

For classical systems, this entropy takes the form 

Sbg = —k J dxp{x) h\p(x) ( J dxp(x) = 1), (3) 

where x denotes a generic point in the full phase space of the system (Gibbs' T space, having typically a 
dimension equal to 2dN, where d is the space dimension where the system evolves). This expression cannot 
be used for too thin distributions, since they would lead to negative entropies. For example, if p{x) = 1/a 
(a > 0) within a one-dimensional (dimensionless) interval whose width is a, and zero otherwise, the entropy 
will be given by Sbg = A; In a. Clearly, this expression is thermodynamically non admissible for a < 1. This 
is a well known fact, and is currently considered to reflect the deep quantum nature of the real world. 
For quantum systems, the Boltzmann-Gibbs entropy takes the form 



Sbg = —kTrplnp (Trp = 1) 



(4) 



where p denotes the density matrix of the system. In this form, the BG entropy is frequently referred to 
as the von Neumann entropy. 

The BG entropy and its associated statistical mechanics provide an extremely useful tool fo r studying a 
wide variety of physical syste ms. However, not all (see, for instance, Gibbs . 19021 : Fermi . 19361 ]). Indeed, as 
suggested since at least 1988 [Tsallisl . 19881 ] . more general or different entropic functionals become necessary 
for the statistical mechanics of other, more complex, systems. 
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The systems for which it is certainly appropriate to apply the BG entropy and consistently associ- 
ated theories can be looseljQ characterized by short-range space-time correlations, Markovian processes 
(short memory), additive noise, strong chaos (positive maximal Lyapunov exponent), ergodic dynamics, 
continuous (Euclidean or Riemannian) geometry for the dynamical occupation of phase space, short-range 
many-body interactions, weakly quantum-entangled subsystems, linear/homogeneous Fokker-Planck equa- 
tions, Gausssian distributions. Such systems neatly benefit from the additivit^ of the BG entropy, and 
typically yield exponential dependences. For example, to start with, the maximization of Sbg under ap- 
propriate constraints yields, for the canonical ensemble, the celebrated BG weight itself: 

p (BG) K e - m t (5) 

where Ei is the energy of state i of a Hamiltonian system satisfying specific boundary conditions, and 
3 = 1/kT. 

There are, however, complex natural, artificial and social systems which, in contrast with the above, 
can be loosely characterized by long-range space-time correlations, non-Markovian processes (long mem- 
ory), additive and multiplicative noises, weak chaos (vanishing maximal Lyapunov exponent), nonergodic 
dynamics, hierarchical (typically multifractal) geometry for the dynamical occupation of phase space, long- 
range many-body interactions, strongly quantum-entangled subsystems, nonlinear/inhomogeneous Fokker- 
Planck equations, non-Gausss ian distributi ons. A quite wide class among these (though surely not all) can 
be handled with the entropy (Tsallisl . Il988| 

l-T w v q W 

S q = k ^f^- (22p i = l;q€n;S 1 =S BG ) i (6) 

^ i=l 

and its analogous continuous and quantum versions. For the particular case of equal probabilities we obtain 

S q = kln q W, (7) 



where 



z x ~ q - 1 

ln q z = — (z G 7£_t_;lni z = Inz) . (8) 



By using this function, entropy ([6]) can be conveniently rewritten as follows: 

W 1 

Sq = k y~] Pi hl q — (9) 

If A and B are two probabilistically independent systems, we straightforwardly verify that 
S q (A + B) = Sq(A) + Sg{Bl _ 5^4) S^) ^ 

If If lp la 

which is a direct consequence of the property 

ln q (xy) =ln q x + ln q y + (1 - q)QiL q x)(]n q y). (11) 

Therefore S q is nonadditive for q ^ 1 (expressions such as subadditive and superadditive are occasionally 
used to refer to the q > 1 and q < 1 cases). 

The respective domains of applicability of Sbg and of S q (q ^ 1) can be simply characterized through 
the equal probability case of a large system, i.e., having a number of elements N » 1. More precisely, if 
W(N) oc fj, N (fi > 1), thermodynamical extensivity (i.e., S(N) tx N in the N — > oo limit) is satisfied by Sbg- 
If, in contrast, we have a system constituted by strongly correlated elements such that W(N) cx N p {p > 0), 
we verify that S q (N) oc N for q = 1 — 1/p < 1. In other words, it is the requirement that the entropy 



loosely we refer to the fact that, in amazingly many cases, necessary and/or sufficient conditions are not available on 
rigorous grounds. 

2 An entropy 5 is said additive [Penrose! . Il97dl ] if, for two probabilistically independent systems A and B (i.e., for discrete 
systems, pf+ B = pfpf , V(i, j)), S(A + B) = S(A) + S(B). This property is straightforwardly satisfied by S BG . 
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satisfies thermodynamic extensivity which determines the appropriate value of the index q to be used. The 
main properties of such systems typically exhibit asymptotic power-laws. For example, in what concerns the 
distribution of energies of a canonical system (i.e., in the rmal cont act with a thermos t at), t he extremization 
of S q under appropriate constraints [Tsallis et al . 19981 ] (see also (Curado Tsallisl 199 yields 



P 



(?) 



ex e, 



-PEi 



(P 



(1) 



oc e 



(12) 



where (3 is related with the temperature T, and the (/-exponential function is defined as the inverse of the 
q- logarithmic function 0, i.e., 



e z q = [1 + (1 - g)zj_ 

with [(...)]+ = (...) if (...) > 0, and zero otherwise. 
Let us notice that, if we extremize say S2 



1/(1-9) 



(z eTZ;el = e 2 



q instead of S q , we obtain 



P 



(?) 



oc e 



2-q 



(13) 



(14) 



which interchanges the q < 1 and the q > 1 intervals. We come back onto this point la ter on. 



The nonadditiy e entrop y S q and its associated nonextensive statistical mechanics [Tsallis Sz Brigattil . 



op > 

20041 : iTsallid . l2009al E | 2010( | have already received a large number of theoretical, experimental, observational 



and computational verifications, including the prediction of the index (indices) q from first principled for 
some systems, and its prediction in terms of mesoscopic quantities for some others. Nevertheless, several 
interesting, delicate and/or elusive points are still open to basic research. It is the purpose of the present 
paper to review several (inter-related) among them, and to indicate their present status of understanding. 

Before introducing the present list of open points, and in order to clarify what kind of natural, artificial 
and social phenomena we have in mind, let us mention several of such phenomena that have already been 
reported in the literature in the frame of the nonadditive entropy S q and its associated nonextensive 
statistical mechanics. 

Among others we have (i) The velo city distribution of (cel ls of) Hydra viridissima follows a q = 3/2 
probability distribution function (PDF) [Upadhvava et al . 200l[ |; (ii) The velocity distribution of (cells of) 
Dicty ostelium d i scoide um follows a q = 5/3 PDF in the vegetative state a nd a q = 2 PDF in the starved 



state [Reynolds . 2010j ]; (iii) The velocity distribution in defect turbulence Daniels et all . 2004]; (iv) The 
velocity distribution of cold at oms in a dissipative o ptical lattice [Douglas et all . 120061 ] ; (v) The velocity dis- 
tribution during silo drainage (Arevalo et al . 2007a b|; (vi) The velocity distribution in a driven-dissipative 
2D dusty pl asma, with q = 1.0 8 ± 0.01 and q = 1.05 ± 0.01 at temperatures of 30000 if and 61000 K 
respectively Liu fc Gored . 2008]; (vii) The spa tial (Monte C arlo) distributions of a trapped 136 Ba + ion 
cooled by various classical buffer gases at 300 K [DeVod . l20Qg| |: (viii) The distribution s of price returns and 
stock vo l umes at the stock exchange, as well as the volatility smile [Borland . l2002al lbl: [5 sorio et all . |2004j; 
Queirosl . 20051 ]; (ix) The distributio ns of returns of magnet ic field fluctuations in the solar wind plasm a 



as observed in data from Voyager 1 Burlaga &; Vinas . 2005 ] and from Voyager 2 |Burlaga &; Nesd . 20091 ]: 
(x) The distributions of returns in the Ehrenfest's dog-flea model iBakar & Tirnakll 120091 . l20ld |: (xi)The 
distributions of returns in the coherent noise model Celikoglu et al . 2010l ]; (xii) The distributions of returns 
of the a valanche sizes in th e self-organized critical Olami-Feder-Christensen model, as well as in real earth- 
quak es (Caruso et al 20071 ]: (xiii) The distributions of angles in the HMF m odel Movano fc Anteneodol . 
2006| ]: (xiv) The distribution of stellar rotational velocities in the Pleiades Carvalho et al . 20081 ]: (xv) 



The relaxation in various paradigmatic spin-glass substances through neutron spin echo experiments 
Pickup et all . 20091 ]: (xvi) Various properties directly related with the time dependence of the width of 



the ozone layer around the Earth Ferri et al . 2010j ] : (xvii) The distribution of transverse momenta in 
high energy collisions of electron-positron, proton-proton, and heavy n u clei ( e .g., Pb - Pb a n d Au-Au) 
Bediaga et al 1200(1 Iwflk fc Wlodarczvkl . 120091 : iBiro et al 120091 : ICMSll . I2OI0I : ICMSl 120101 : IPHENIXl . 



I j |l> 1 'II 1 1 I u , tl 1 pMMii aiiiia" JBHi^ , , 7 I 17 1 ij ■ — 1 ? 

2010J; IShao et all . |2010j], the flux of solar neutrinos Kaniadakid . Il996l ] . and the energy distribution of 



By first principles we mean from the set of probabilities of the microscopic configurations and its corresponding dynamics. 
For example, whenever the system is a mechanical conservative one, from first principles it is meant from the Hamiltonian or 
from the Lagrangian of its elementary constituents. 
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namical systems Lvra &: Ts allis 



cosmic rays Tsallis et all. I2003H ; (xviii) Various properties for conservative and dissipative nonlinear dy- 



19981: iBorges et all 120021: lAnanos fc Tsallisl. 120041: iBaldovin fc Robledol . 



20041 : iMavoral Robledo 1 . 120051 : IPluchino et alTEoO T. 2008: iMiritello et all. [20091: ; iLeo et at l20ldl : (: 



xix 



- tj , u , _ _ - — i; i m j ■ =^ — ssaj i ij i U ■ 1 __^ iri64 , , — |^^^^^^ ^ < / 

The degree distribution of (asy mptotically) scale-free networks [White et all . 120061 : iThurner et al l2007i | ; 
(xx) Tissue radiation response Sotolongo et al 2O10l ] : (xxi) Overdamped motion of interacting particles 
Andrade et al . l2010l pl. 



2. Some Open Points 

2.1. q-generalized Lyapunov spectrum and Pesin identity 
2.1.1. q-generalized Lyapunov spectrum 

Let us illustrate this point through a simple example, namely a one-dimensional dissipative map, e.g., the 
z-logistic family 

x t +i = 1 - a\x t \ z (t = 0,1,2,3,...; -1 < x t < 1; z > 1; < a < 2) (15) 

The particular case z = 2 corresponds to the standard logistic map; the z — > 1 limit corresponds to the 
tent map. These maps metrically, but not topologically, differ for different values of z. For the simple z = 2 
case, and with y = x + 1/2, we obtain the traditional form 

y t+l =ny t (l -y t ) (0</i<4;0<y 4 <l). (16) 

For a increasing above zero, a succesion of fixed points and fixed cycles occur, separated by doubling- 
period bifurcations. These bifurcations accumulate as a approaches a special point, a c (z), the first edge of 
chaos. For z = 2 it is o c (2) = 1.40115518909... 

The sensitivity to the initial conditions £ for a one-dimensional dynamical system is defined as follows: 

Ax(t) 



lim 



m = i) 



(17) 



Aa;(0)^0 Ax(0) 

where x denotes the phase space variable. The system is said strongly chaotic (or simply chaotic) if £ 
exponentially diverges with time. In such cases we can define the Lyapunov exponent A through 

(18) 



or, more precisely, through 



A 



£(i) ~ er 

ln£(i) 



lim 

t— ¥OC 



(19) 



At the edge of chaos, A vanishes and £ increases slowly with t, in fact algebraically at large enough 
values of t. For a increasing above a c (z), A greatly oscillates in a complex manner, being however positive 
for most of the values of a. The sensitivity £ is, in fact, quite generically expected to satisfy 

<K K*J qsm (£(0) = 1), (20) 



dt 



hence [Tsallis et al Il997l : IBaldovin and Robledol . l2002al lbl: iRobledol . [2006] 



-'Qsen 



+ 



(21) 



where q sen = 1 if the Lyapunov exponent Ai = A 7^ (strongly sensitive if Ai > 0, and strongly insensitive 
if Ai < 0), and q sen ^ 1 otherwise; sen stands for sensitivity. At the edge of chaos, q sen < 1 (weakly 
sensitive), and at both the period-doubling and tangent bifurcations, q sen > 1 (weakly insensitive). The 
case q sen < 1 yields, in (|2"Tj) . a power-law behavior £ oc tV( 1 ~ 9se '0 j n the limit t — > 00 . This power-law 
asymptotics were since long known in the literature Grassberger and Scheunertl . Il98ll : ISchneider et al . 



4 Many other ph enomena have been looked a t along similar lines (e.g., biological evolution [Tamarit~et al, 1998], turbulence in 
electron plasma [Anteneodo fc Tsallisl . Il997j |). 
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19871 : lAnania fc Politl Il988l : lHata et all . Il989l : iMori et all . Il989f ] . The case Qsen *C 1 is in fact more complex 
than indicated in Eqs. (|20j) and (j2~Tj) . These equ ations only reflect the maximal values of an entire family, 
fully (and not only asymptotically) described in jRobledol . 2006 ; Mavoral &: Robledo . 20051 ]. 

The rigorous necessary and sufficient conditions for behaviors such as those indicated in Eqs. \2(fy and 
(21)) . for generic conservative and dissipative nonlinear maps, stand at present as an open problem. 



2.1.2. Pesin-like identity 

Let us now focus on the time evolution of the entropy of the above illustrative maps. We shall illustrate 
with the map defined in Eq. (|15p . We divide the phase space x E [—1, 1] in W equally spaced little intervals 
denoted with i = 1,2,3, W^(with W >> 1), choose one of them, and put within M (randomly chosen) 
initial conditions (with M » 1; typically M ~ 10 x W). At time t we have in the i th interval Mi(t) 
points; naturally Yli=i -^i(^) = M. We then define a set of probabilities through Pi(t) = Mi(t)/M (hence 
^j = i Pi (t) = 1), which enables the calculation of the entropy (we take k = 1) 



i 



i 



We next define the entropy production (per unit time) as follows: 

S (t) 

K Q = lim lim lim — - — 

t-s>oo W^foo M-s-oo t 



(22) 



(23) 



There typically exists an unique value of q, noted q pro duction, such that K qprodnotion is finite. For q > qproduction 
(q < qproduction), Kq vanishes (diverges). For various one-dimensional dissipative maps (and also two- 
dimensional conservative ones) we verify that 



and 



Qproduction Qsen 



^^Qproduction ^Qscn 



For strongly chaotic maps, i.e. with Ai > 0, we verify that 

Qproduction Qsen 1 i 



and 



K 1 = X 1 > 0. 



This last equality can be seen as a Pesin-like one. 

For weakly chaotic ones (typically at the edge of chaos), i.e., with Ai = 0, we verify that 



and 



Qproduction Qsen *^ 1 j 



^ Qproduction ^Qsen ^ ^' 



(24) 
(25) 
(26) 
(27) 

(28) 
(29) 



This last equality can be seen as the ^-generalization of the Pesin-like equality. See Fig. [2] 

The rigorous necessary and sufficient conditions for equalities such as those indicated in Eqs. \2J$ and 
A25\), for generic conservative and dissipative nonlinear maps, stand at present as an open problem. 



2.2. Geometry of occupation of phase space, Hilbert space 
2.2.1. Additivity versus extensivity of the entropy 

The additivity of an entropy only depends on the mathematical functional which expresses the macroscopic 
(thermodynamical) entropy in terms of its basic probabilities. Therefore, as already said, the BG entropy 
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TIME 
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• a=1.6 
a a=2 

- fit slope 0.36 
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Q. 
O 
DC 



W 
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20.0 30.0 
TIME 



40.0 



Fig. 2. Time evolution of the BG entropy for typical values of a of the z — 2 logistic map, corresponding to positive Lyapunov 
exponents. Top: We notice that only for q — 1 we observe a linear intermediate behavior of the entropy before saturation. 
This linear region is larger for larger W . Bottom: The entropy p roduction Krq pe r unit time decreases when the value of a 
corresponds to a smaller value of the Lyapunov exponent. From [Latora et all . 120001 ] . 



Sbg is additive. So is the so-called Renyi entropy (useful in the geometrical description of hierarchical 
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structures such as multifractals) , denned (for the discrete case) as follows 

S*= \: q J (qen;S« = S BG ). (30) 

Indeed, it can be straightforwardly proved that, if A and B are two probabilistically independent systems, 
then S*(A + B) = Sf{A) + S*{B). 

In contrast, the extensivity of an entropy is a more subtle concept, and it depends not only on its 
functional form but also on the specific system under consideration, i.e., on its microscopic probabilistic 
correlations. Consequently, while the expression "the BG entropy is additive" is definitively correct, the 
expression "the BG entropy is extensive" is incorrect. The correct expression would be "the BG entropy 
of this class of systems is extensive" . 

Consider a system S = A\ + A2 + ... + An constituted by N (not necessarily independent) identical 
elements or subsystems {Aj}. An entropy S of that system is extensive if < limjv^oo S ffl < 00, i.e., if 

S{N) ocN (N -»• 00) . (31) 

The important difference between additivity and extensivity can be illustrated through probabilistic 
systems of N identical binary random variables. If the variables are distinguishable, we may represent them 
in the following triangular form: 

(N = 0) lxl 

(N = 1) 1 x r\fi 1 x ri ; i 

(N = 2) 1 x r 2j0 2 x r 2 ,i 1 x r 2i2 

(N = 3) 1 x r 3j0 3 x r 3j i 3 x r 3i2 1 x r 3)3 

(N = 4) 1 x r 4i0 4 x r 4 ,i 6 x r 4)2 4 x r 4)3 1 x r 4j4 



where < r^^n < 1 and 



n=0 v ; 



If the distinguishable variables are independent we obtain the particular case 

r N ,n=p N - n (l-p) n (0<p<l) (33) 



(AT = 0) lxl 

(AT = 1) lxp lx(l-p) 

(AT = 2) lxp 2 2xjj(l-p) lx(l-p) 2 

(AT = 3) lxp 3 3xp 2 (l-p) 3xp(l-p) 2 lx(l-p) 3 

(AT = 4) lxp 4 4xp 3 (l-p) 6xp 2 (l-p) 2 4xp(l-p) 3 lx(l-p) 4 



If p = 1/2, this triangle becomes 



(N = 0) lxl 

(AT = 1) 1 x \ 1 x \ 

(N = 2) lxl 2x1 i x i 

(N = 3) lxl 3 X i 3 X i i x i 

(N = 4) Ixl4xi6xi4xllxl 
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(N = 0) lxl 

In = i) i x ± i x i 

(JV = 2) lxi 2xi lxi 

(N = 3) lxi 3xi 3xi lxi 

(TV = 4) Ixi4xi6xl4xilxi 



In the presence of correlations, the generic triangle can take various forms characterized by the sets 
{ r N,n} i ViV. For example we have the Leibnitz triangle. Its generic term is given by 

rN -= m N + i - (34) 

The Leibnitz triangle satisfies the following remarkable property (from now on referred to as the 
triangle Leibnitz rule): 

rN,n + rN,n=l = TN-l,n ■ (35) 

This property implies probabilistic scale-invariance. Indeed, it implies that the marginal probabilities of a 
iV-system coincide with the joint probabilities of a (N — l)-system. 

Let us now turn back to the entropy. For a wide class of systems (which includes all the above prob- 
abilistic triangles) a value of q exists, noted q en tropy, such that S qentropy is extensive. For all the above 

examples it is Qentropy = 

1, i.e., Sbg(N) oc N (N -»■ oo). 
Let us now severely restrict the admissible probabilistic region (defined as the set of configurations 
whose probability is strictly positive). We are referring to cases where the total number, noted W e ff(N) 
(eff stands for effective), of admissible configurations is much smaller than W{N) for TV >> 1. In other 
words, cas es where lim y y_ > . 00 W P ff(N)/W(N) = 0. Two such examples are indicated in what follows (see 
details in (Tsallis et all . [2005] ) . These triangles have nonzero probabilities only along a (left) strip whose 
width is denoted d+1. The two examples here respectively correspond to d = 1 and d = 2. 

1 x 1 
1 x | 1 x \ 
1 x \ 2 x ± 1x0 
1 x ± 3 x ± 3x0 1x0 



(N = 


0) 


(N = 


1) 


(N = 


2) 


(N = 


3) 


(N = 


4) 


(N = 


0) 


(N = 


1) 


(N = 


2) 


(N = 


3) 


(N = 


4) 



2 

1 x | 4 x I 6x0 4x0 1x0 



1 x 1 

1 x \ 1 x \ 

1 X i 2 X i 1 X | 

lx|3x^ 3x| 1x0 
lx I 4x I 6x I 4x0 1x0 



It can be shown that, for this family of triangles, we have W e jf(N) << W(N) = 2 (N » 1), and 
consistently 

Qentropy = 1 — ^ (d = 1, 2, 3, ...) . (36) 

2.2.2. Some many-body physical examples enabling first-principle calculations of q 

Quantum entanglement is caused by the intrinsic nonlocality of quantum mechanics. This nonlocality makes 
the elements of a iV-body (iV >> 1) system to be strongly correlated, which diminishes considerably the 
size of the admissible space of microscopic configurations. Depending on the specific system, Sbg(N) 
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might be extensive. We consider here what is currently referred to as the block entropy, i.e., the entropy 
of a subsystem which is very large and nevertheless much smaller that the entire iV-system. The elements 
of this subsystem are quantum-entangled with all the ./V elements, but our detector only detects (for 
whatever reason) the elements of the subsystem. For the subsystem, it frequently happens that its Sbg 
is nonextensive, and therefore inadequate for thermodynamical purposes. We review here two magnetic 
examples for which a nontrivial value q en tropy exists such that the block entropy S qentropy is extensive, thus 
reconciling the subsystem entropy with thermodynamics. 

Both examples concern chains of many-body spin systems with short-range interactions at T = 0, i.e., 
at their fundamental state. The total iV-system is in a pure state, and therefore its entropy vanishes. A 
block of successive L spins (with 1 << L « N) is, however, in a mixed state and therefore its entropy 
S q (L) is different from zero. By S q (L) we mean precisely 



S q {L)= l -^f, (37) 



where 

PL = Tr N _ L p N , (38) 

Pn being the density matrix of the iV-system. We have that Trp 2 N = 1 (i.e., a pure state, hence S q {N) = 0), 
whereas Trp\ < 1 (i.e., a mixed state, hence S q (L) > 0). 

The first example is a pure ferromagnet with spin 1/2 anisotropic first-neighbor interactions in the XY 
plane in the presence of a transverse magnetic field at its critical value (a T = quantum phase transition 
from the paramagnetic to the ferromagnetic state). It contains the I sing (XY isotropic) ferro magnet as the 
central charge c = 1/2 (c = 1) particular case. We analytically find Caruso Tsallisl 20081 ] that 



V9 + c 2 - 3 

Qentropy = ~ (c > 0) . (39) 

The second example is a random magnet of the Heisenber g type with spins S, in the absence of any 
external field. It is numerically found Saguia &; Sarandv . 2010l | that 



1.67 1.67 . 

9^ = 1- — = 1-^2^1) (C " 0) - (40) 



Both expressions for q en tropy ar e represented in Fig. [3J 




Fig. 3. Curves of q versus 1/c. The red dot corresp onds to the BG limit. T h e continuous (dashed) cur ve corresponds to the 
T = linear chain for the pure (random) magnet in [Caruso ~H Tsallisl . l200j ] ( [Saguia fc Saran dvl. l2010l ]). 
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Since the block is not in a pure but in a mixed state, it necessarily has a nontrivial energy distribution. 
Which distribution? It could be given, in the N — > oo limit, by say 

PL^e-±Z e y r3yHL (L»l), (41) 
where q en ergy could well be a (relatively simple) function of q en tropy- 

Qenergy — ^{.Qentropy) • (^) 

In the c — > oo limit one naturally expects the BG limit q en ergy = qentropy = 1 to hold (i.e., h(l) = 1). 

The block density matrix pi stands at present as an open problem. If it turns out to be of the form 
l4l\) , what would be the values of q en ergy and of (3 qenergy ? 



2.3. q- generalized central limit theorems 

2.3.1. q-product 

For reasons that wil l soon become clear let us introduce a ^-generalization of the product, the q-product 



Nivanen et al . 120031 : iBoreesl . 12003 ] : 



x® q y = [x 1 q + y 



,1-91- 



i (x ®i y = xy) . 



(43) 



The main reason for defining such a product is that it satisfies the following remarkable property: 

hx q (x® q y) = ln q x + ln q y. (44) 

Similarly to Eq. (Illh . which essentially reflects the nonadditivity of the entropy S q , Eq. (I44p reflects its 
possible extensivity in the presence of strong correlations characterized by q. Let us illustrate this fact. 
Assume that we have a system whose subsystem A (B) has N A >> 1 (N B >> 1) strongly correlated 
elements such that its nonzero-probability configurations are equally probable, and equal to 1/W(N A ) 
(1/W(N B )). Assume also that the correlations are such that W(N A ) oc (N A ) P (W(N B ) oc (N B ) P ) with 
p £ 1Z. Then we straightforwardly verify entropic extensivity for q = 1 — 1/p, i.e., that 5 1 _ 1 / p (Aa + N B ) oc 
N A + N B . 

We shall focus here on the (/-product for the case q > 1, and x and y non-negative real numbers (see 
Tsallis fc Queirosl . l2007l | for further details). It is straightforward to check that this product is commutative, 



associative, it has an inverse, a unit, and a zero. But there are preliminary indications that there is 
no associative generalized sum with regard to which the g-product is distributiv e. The argume nt goes 
essentially as follows. We first develop the g-product as follows (see Appendix A in [Tsallisl . 12009^ : 



x® q y = xy^l + (q- l)(lnx)(lny) + ^(q - l) 2 [(In 2 x)(lny) + (lnx)(ln 2 y) + (In 2 x)(ln 2 y)] + ... j . 
We next assume the existence of a generalized sum x® q y whose development would be as follows 0: 
x(B q y = (x + y)|l + {q - l)ai(x,y) + ^(q - l) 2 a 2 (x,y) + ...j , 



(45) 



(46) 



where ai(x,y) = a\(y,x) and a2(x,y) = a2(y,x) are functions to be found by simultaneously imposing (i) 
distributivity of the q-product with regard to x © g y, and (ii) associativity of this generalized sum. In other 
words we impose 



and 



X <S>g (y ® q Z) = (X (g)g y) ©g (X (g)g z) 



x(B q (y® q z) = {x(B q y)e q z. 



(47) 



(48) 



We use x ® q y instead of simply x @ q y because, for different purposes, the latter is already defined in the literature. It is 
called g-sum and is defined as x ® q y = x + y + (l — q)xy. Let us say right away that x ® q y is not distributive with regard to 
x ® q y, as can be easily verified. 
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Also, we impose x ® q = 0, V</, which immediately implies a±(x, 0) = a,2(x, 0) = 0. By developing in powers 
of (q — 1) both relations (|47l) and (|48l) . we obtain strong cons traints on a-\ ( x, ?/) and d2(x,y). Although 
the full discussion remains to be done, a preliminary analysis [Curado et all . 120081 ] suggests that no such 
ai(x,y) and a2{x,y) can exist. If indeed they do not exist, no algebra exists whose product would be the 
g-product. That would be some kind of new mathematical structure, satisfying all the axioms of an algebra 
but distributivity. 

The distributivity of the q-product stands at present as an open problem. Depending on whether it exists 
or not, and on what would be the analytical definition of x® q y, is it possible to construct some structure 
having some similarity with a vector space? 



2.3.2. q-Fourier transform and discussion of its inverse 

By using the definition (143p we can define the following (/-generalized Fourier transform (q-FT) 
Umarov et al 120081 . l2010l ]: 



F q [f(x)}(0 = I dxe** ® f(x) (q > 1) (49) 
For non-negative real f{x) we have that 

F g [mm = I dxe^VW-'fix) (q > 1) (50) 



It is clear that 1-FT recovers the standard Fourier transform (FT), which is a linear integral transform. It 
is also clear that, for q ^ 1, the q-FT is a nonlinear integral transform H. It has been so defined in order 
to be closed with regard to the family of the (/-Gaussian distributions. Let us be explicit. We consider the 
following g-Gaussian 



f(x) = G q (x) ^ - [1 ■ (1 q)Px ] \ q _ , = „ 6g _ py2 ! ? > (»: c, < 3) . (51) 



Jdy[l-(l-q)/3y2}— fdye q 



where we remind that the support is infinite for 1 < q < 3, and is finite for q < 1 □. Eq. (15ip can be written 
explicitly as follows: 

G q (x) =N qV ^[l-(l-q)(3x 2 }^ (p>0;q<3), (52) 



6 If i is a variable which carries physical dimensions, it is convenient to define definition (|49[) as follows: F q [f(x)](£) = 

fd(xf )ef 6 ® (f(x)/f ) = fd(xf Q )e l f [f( - x)/fo]q '\f(x)/fo) = f dxjf [1(x)/ ^ f{x) (q > 1), where f > is a 
reference value. If /(0) is finite, a simple choice would be to just adopt fo = /(0). With this generalized definition it follows 
that x£ is a pure number, which physically is very convenient. If x is already a pure number, we can of course adopt /o = 1. 
7 g-Gaussians are referred to with various names in the literature. In the form f(x) oc l/[a 2 +x 2 ] K , they are occasionally called 
"generalized Lorentzians" or "kappa distributions" in areas such as plasma physics, and "Barenblatt form" in the area of 
porous media. Also they recover the Student' s t-distributions and the r-distributions for special values of q. Because of that, 
in finance, they are loosely referred to as "Student' s t-distributions" even outside these special values of q. 
8 Let us remind that the g-Gaussian form is normalizable only for q < 3. Its variance is finite for q < 5/3, and diverges for 
5/3 < q < 3. Its q- variance, however, remains finite for any q < 3, i.e., as long as it is normalizable. 
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with 



N q = ( 



q-1 



7T 



1/2 



q-1 



TT 



3-q 



7T 



3-g 
2((?-l) 



3-g 
1/2* \ 2(l-q) 



1-5 



if 1< g < 3 . 



if g=l, 



if g < 1 , 



(53) 



We can verify that, for q > 1, 



with 



and 



F q [G q (x)](0 = e/ e , 



1 + 



* 3-q' 
(3 - q)N^ 



This expression can be conveniently rewritten as follows: 

0)^^- r(3-^ 2(1 - s) 



(54) 
(55) 

(56) 
(57) 



For q = 1, we recover /3/3 = 1/4, which leads, in quantum mechanics, to the Heisenberg uncertainty 
principle. Notice also that Eqs. (|55|) and (|56|) are invertible, i.e., (q,/3) — > (q,/3) and (q,(3) — > (g,/3). 

The above relations show that the q-FT is invert ible with i n the class of the g-Gaussians. It is not 
so in general. Indeed, it has been shown by Hilhorst Hilhorstl . I2009I . l201d( | that, for a given value of q, 
one-parameter families of functions {f(x)} exist such that their q-FT does not depend on that parameter. 
Therefore the q-FT has not always an unique pre-image, i.e., although the g-FT is invertible within the 
closed class of g-Gaussians, it is not invertible in general. This peculiar property does not exist for q = 1, 
but it does emerge for q > 1. Therefore an interesting question can be put: Does a procedure exist which 
given a (say non-negative, for simplicity) function f(x), enables us to go forward an d backward through 
the q- FT and come back to the same function? The answer is yes, as recently shown in [jauregui <fe Tsallis . 
2010a]] . The procedure is as follows. 

Given a non-negative function f(x), we make the transformation x — > x + y, and then calculate the 
q-FT on the variable x, i.e., we calculate 



F q [fyU,y) 



+oo 



f(x + y)ef x ^ x+ y^ q ~ 1 dx 



(58) 



Next, by using the recently introduced g-Dirac delta [Jauregui fc Tsallisl . [2010a] ] . we can recover the function 
f(y), or equivalently f(x), through 



2-q 
2vr 



2-q 



(59) 



for all points of the interior of the support of f(x) (see [jauregui &; Tsallis] . 201 Obi ] for details concerning 
the points at the edge of the support). 
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2.3.3. q-independence 

Two random variables X and Y, with respective distributions fx(x) and /y(y), are said independent if 
the joint distribution satisfies 

f(x,y) = f x (x)f Y (y) (60) 
This implies that the distribution of their sum is given by 



fx+v{z) = J dx J dyf(x,y)6(z - x - y) = J dxf x (x)f Y (z - x) , (61) 
where 5 denotes the Dirac delta. If we take the FT we obtain that 

F\f x+ Y(*M) =F[f x (z)](0F[fY(z)](0- (62) 
Due to this property, we can alternatively define independence of X and Y by saying that the FT of 
/x+y(z) equals the product of the FT of fx(x) and the FT of /y(y). If the variables are not independent 
(i.e., if they are correlated), we have that 



F[fx + Y(z)](t) = j dx J dyf(x,y)5(z-x-y)^F[f x (z)](0F[f Y (z)}(0. (63) 
We shall next define a special type of correlation between X and Y, named q-independence, by imposing 

F q [f x+Y (z)](0 = F q [f x (z)W) ® f F,LM*)](£) . (64) 

with q given by Eq. ()55j) . This definition might seem rather esoteric at this stage, but it will turn out to be 
amazingly frequent in nature, as we shall see later on. Of course 1-independence recovers independence. 

2.3.4. q- generalized central limit theorems 

In the empirical sciences, the repetition of an experiment (in nearly equal conditions) is a must if we wish 
to increase the experimental precision. If we do it ./V times, we then consider as valid experimental value 
the arithmetic mean of those N results. This ubiquitous fact points towards the importance of considering 
random variables of the type Sn = X\ + X2 + ... + Xn, which constitutes the scope of the Central Limit 
Theorem (CLT). This important theorem admits a variety of forms, but basically it states that, if {Xi} 
are equally distributed independent (or quasi- independent in some sense), have a finite variance, and N is 
increasingly large, the sum Sn always converges, after appropriate centering and scaling, onto a Gaussian, 
which is therefore called the attractor of the sum. If the variance diverges, the attractor becomes instead 
a Levy distribution, also called a-stable distribution with a < 2 (this is sometimes referred to as the 
Levy-Gnedenko CLT). 

It is quite natural that the attractors are neither Gaussian nor Levy distributions if strong correlations 
are present between the {X{\ variables. It has been recently shown Umarov et all . 2008] that, if the variables 



are equally distributed q-independent and a conveniently (/-generalized variance cjq is finite, the N — > 
00 attractors are (/-Gaussians. If that (/-generalized variance diverges, then th e nature of the attr actors 
is different and they are referred to as the (q, a)-stable distributions L q ^ a (x) (Umarov et all . I2OI0I ] . The 
ct — )■ 2 limit of the {q, a)-stable distributions are the (/-Gaussians L q ^{x) = G q (x); the q — > 1 limit of 
the (q, a)-stable distributions are the a-stable distributions L\^ a {x) = L a (x); the limit L\^{x) = G{x) 

f°° dxx 2 [f(x)] Q 

corresponds to Gaussians. Q is defined as Q = 2q — 1 (with 1 < q < 2), and ctq = ySL dx \j^yp ■ These 
theorems are schematica l ly depicted in the follow ing Table (where the C q> aS are positive coefficients). See 



Tsallis fc Queirosl . 120071 : iQueiros fc Tsallisl . 120071 ] for typical illustrations of the four types of attractors. 

There are many Gaussians in nature: it seems reasonable to believe that this is due to the CLT, which 
shows that details can be of no importance if the number of involved random variables is large. Similarly, 
and for the same reason, we also expect many (/-Gaussians to exist in nature (as well as in artificial 
and social systems), as long as ^-independence emerges naturally in many complex systems, where strong 
correlations between elements is an important ingredient. And it is precisely this situation that has been 
profusely illustrated in Section 1. 
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g = 1 [independent] 


g / 1 (i.e., Q / 1) [globally correlated] 


<7Q < OO 

(a = 2) 


G(x) 

[with same <Ti of /(a;)] 


G 9 (x) 

[with same <tq of /(x)] 
G g (x) ~ G(x) if x << x c (q, 2) 
G,(x) ~ C q>2 /\x\ 2 /(«-V if | x | » Xc (q,2) 
for g > 1, with lim g _^.i x c (q, 2) = oo 


C7Q -> OO 

(a < 2) 


L a (x) 

[with same x — > oo behavior of f(x)] 

L a {x) ~ G(x) if \x\ « x c (l,a) 
L a (x) ~ Ci )Q / x 1+a if x >> x c (1,q) 
with lim a _ 5>2 ^c(l) a ) = 00 


[with same x — > oo behavior of /(x)] 

i g ,a~^:2 terme * ate) /ki ( q) 

if i^g, a) << x << Xc 2 ^, a) 
L q>a ~c(% stant) /\x\^^ 

Til (2)/ \ 

if x >> Xc ((/, a) 



Let us incidentally mention that the lack of inverse of the q-FT has made Hilhorst [H ilhorstJ . [2009, 

20101] to disregard q-Gaussia ns as attractors. A detailed reply to his claim has been recently made available 



m 



Umarov fc Tsallial . l20ld ]. which reinforces that g-Gaussians are, in many respects, very special distri 



butions. It must be also taken into account that, from a different perspective, a theorem like the g-CLT 
(and even other forms as s ociat e d with other type s of correlations) has been also proved without using the 
q-FT jVignat fc Plastinol . 120071 : lHahn et all . l2010l ] . 

The above discussion main ly focused q > 1. However, the q-FT has also been addressed for q < 1 (see 



Nelson fc Umarovl . 120081 . l2010l ]). 



As an interesting mathematical challenge we may now ask: Is it possible to alternatively prove the 
q-CLT by simultaneously using the q-FT and relation i59\) ? 



2.3.5. Possible relation between q-independence and scale-invariance 

Two probabilistic models, [Moyano et a I l2006l ] and IThistleton et a respectively, involving N equally 

distributed random variables were introduced some time ago. Their numerical discussion suggested that, 
in the N — > 00 limit, q-Gaussians emerg ed with q < 1, after appropriate centering and scaling. It was 
however proved Hilhorst &; SchehH . 2007 ] quickly after that the limiting distributions of these two models 
are not exactly g-Gaussians, even if numerically they are amazingly close to them. This interesting result 
put forward a relevant question, which we describe now. 

We will consider scale-invariant a iV-particle probabilistic model which satisfies the following property: 

fN-i(xi,x 2 ,...,x N -i) = / dx N f N (xi,x 2 ,...,x N ) , 



(65) 

where /at(xi, x 2 , xjv) is the joint probability distribution associated with the random variables 
(xi, X2, xjv), satisfying J dx\ dx 2 ...dxN /jv( x l> x 2 ; %n) = 1 (ViV). In other words, a iV-particle prob- 
abilistic model is said scale invariant if the marginal probabilities after tracing over any particle the iV- 
particle joint probabilities coincide with the (N — l)-particle joint probabilities. If the system is made by 
N binary random variables, scale invariance is nothing but what is sometimes refer red to as the Leibn itz 
triangle rule (not to be confused with the Leibnitz chain rulel), i.e. (see for instance Tsallis et al . 20051 ] ). 



r N)11 -i + r N)11 = r7v_i jn (A = 2,3,4,.. 



n 



0,1,2,..., A) 



with 



N 

E 

n=0 



A! 

(A - n) \n\ 



r N>n = l (VA) 



(66) 



(67) 
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It is straightforward to prove that probabilistic independence yields scale-invariance. Indeed, if no 
correlation is present, we have that /n(xi,X2, ■■•,%n) = fi(%i) fi(%2) —,fi(%N) with J dx fi(x) = 1, which 
trivially implies scale-invariance. But, even in the presence of correlations, scale-in yariance is possible . 
The celebrated Leibn itz triangle is one such example. Both models introduced in [Movano et all . 2006; 
Thistleton et all . 2009 ] also are nontrivially scale-invariant. But they are not (/-independent. Indeed, if they 
the N — > oo limiting distributions of both models would be (/-Gaussians, and they are not. 



were. 



We immediately conclude that scale invar i ance i s not sufficient fo r (/-independence. It is nevertheless 
compatible with it. Indeed, in Rodriguez et al . 20081 : Hanel et~a3 . 2009 ] we have introduced scale-invariant 
models whose limiting distributions are (/-Gaussians for values of q both above and below q = 1 (which 
corresponds in fact to independence). 

The present scenario is as follows: scale-invariance clearly is not sufficient for (/-independence, but could 
well be necessary. 

So the following question emerges: What are the exact mathematical implications between (strict or 
asymptotic) scale-invariance and (strict or quasi) q-independence? Or, alternatively, is there some other 
generic condition which, added to scale-invariance, makes it necessary and sufficient for q-independence? 



2.4. q- generalized plane waves, Dirac delta, wave equation and nonlinear 
Klein- Gordon, Schroedinger and Dirac equations 

The (/-exponential function has various remarkable properties which straightforwardly generalize those of 
the exponential function. Let us illustrate this through the following ordinary differential equations (ODE). 
We consider the (linear) ODE 

dy . 
— = aiy (oiGK), 

with y(0) = 1. The solution is given by 

y [x) = e aix . 

We may now generalize Eq. (168p into the following (nonlinear) one: 

dy 
dx 

The solution is now given by 



a q y Q ((a q ,q)£K 2 ). 



y(x)=e q q . 

We may finally unify [Tsallis et al Il999l ] both Eqs. dggj) and {TO} into 



-1- = aiy + {a q -ai)y q ((oi, a q , q) G V?) , 



whose solution is given by 



1--* + 



Q e (l-q)aix 



i 
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(68) 
(69) 

(70) 
(71) 

(72) 
(73) 



This solution makes a crossover from the (/-exponential (for small x) to the exponential function (for 
|(1 — q)a\x\ » 1). As an interesting remark, let us mention that, for q = 2, this solution becomes y= 
- — — q — a q g _ — x , which in the appropriate limits (and with the photonic density of states) recovers the Planck 

law for the black-body radiation! (see details in [Ts allisl . l2009al ] ) . 

Let us now turn onto (/-plane waves and related matters, and later on we will address a possible relation 
between them and the crossover we have mentioned here above. 

To start we remin d that the representation of Dirac delta in terms of plane waves has been recently 
generalized as follows (jauregui &; Tsallis . 2010bl ] in terms of (/-plane waves: 

2 - 



S(x) 



2tt 



dk e, 



-ikx 



(!<(/< 2), 



(74) 
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which implies that a wide class of functions f(x) exist such that 

dx 5{x - x )f(x) = /(s ) (1 < q < 2) . (75) 



(76) 



Let us now focus on the standard one-dimensional linear wave equation 

d 2 $(x,t) 1 d 2 $(x,t) 



dx 2 



dt 2 



for which any function of the type §{kx — ut), twice differentiable, is a solution. In particular, one may 
have a q-plane wave, 

$(x, t) = $ exp ? [i(kx - ut)} ; [$ = *(0, 0)] , (77) 

as a solution of the equation above, provided that jNobre et all . [2010] 

u = ck (l<g<3). (78) 

Its generalization to d dimensions is straightforward (see [Nobre et al . 2010( |). Its plane wave solution 
becomes 



<J?(x, t) = $o exp i(k ■ x — ut) 



(79) 



Let us introduce now the following ci-dimensional nonlinear generalization of the Schrodinger equation 
for a free particle of mass m, 

n 2 _„ [$(x,ty 2 - q 



ih 



0_ 

dt 



1 



q 2m 



-V 2 



(80) 



We notice that the scaling of the wave function by <3?o guarantees the correct physical dimensionalities for 
all terms. This scaling becomes irrelevant only for linear equations [e.g., in the particular case q = 1 of Eq. 
(|80p . N ow, if ones u s es the g-plane wave solution Eq. (I79p by simply replacing k — >■ p/h and u — > E/h, one 
verifies Nobre et a i I2010I ] that the form 



&{x,t) = $ ex P(? 
is a solution of the above nonlinear equation, with 



— {p-x — Et) 
n 



E = p 2 /2m , 



(81) 



(82) 



thus preserving the well known energy spectrum of the free particle for all values of q. 

Let us now add a mass term in the traditional wave equation, and propose the following nonlinear 
Klein- Gordon equation in d dimensions, namely 

2(9-1) 



\ 1 d 2 $(x,t) ... . r . , 
V 2 $(x,t) = — tt~t^ — + q ®(x,t) 



m 2 c 2 



<$>(x,t) 







(83) 



• c 2 dt 2 ' * n 2 

One may verify easily Nobre et all . 2010] that the same (/-plane wave used for the nonlinear Schrodinger 
equation is a solution of Eq. (]83}> . preserving for all q the Einstein relation 

E 2 = p 2 c 2 + m 2 c 4 . (84) 



Let us finally consider a nonlinear generalization of the d = 3 Dirac equation [Dirad . 119281 ] , namely 
Nobre et at l2010f ] 



ih d^j' Q + ihc{a ■ V)$(x, t) = Pmc 2 A i - q) (x, t) $(x, t) 



(85) 



where a x ,a y , a z (written in terms of the Pauli spin 2x2 matric es a T , o, n an d o~ z ) and (3 (written in terms 
of the 2x2 identity matrix /) are the standard 4x4 matrices Libofi, 20031 ] . namely 
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a., 



a x 
a x 



and 



(Jy 
Oy 



/ 
-/ 



O:- 



a y 

Oy 



The new, q-dependent, term is given by the 4x4 diagonal matrix A^\x, t) = 8ij[$j(x,t)/aj] 

Qts (A§\x,t) = Sij). The solution of Eq. (Jg 

four-component column matrix 



(86) 
(87) 

9-1 



where {aj} are complex constants (A^\x,t) = 5^). The solution of Eq. ([85]) we focus on is the following 















$ 3 


%t) 
%t) 




a-2 


exp g 


% -{p-x-Et) 


\$ 4 


%t)) 




W 







<$>(x,t) = 



Substituting this four-component vector into Eq. (|85p . we get, for the coefficients {aj}, precisely the same 
set of four a lgebraic equations corresponding to the linear case, namely (see page 803, Eq. (15.45b) of 



Libofj . l200aj ]). 



(E 
(E 
(E 
(E 



cp z a 3 - c(p x - ip y )a 4 = 

c(p x + ip y )a 3 + cp z a4 = 

mc~)a 3 - cp z a\ - c(p x - ip y )a 2 = 
2^ 



mc )a\ 
mc 2 )a2 

2 



mc 



)a 4 - c(p x + ip y )ai + cp z a 2 = 



(89) 



These equations have, for all q, a nontrivial solution only if the Einstein energy- momentum relation 
Eq. (JHU) is satisfied. 

The above nonlinear Schroedinger equation (I80p can be further generalized as follows: 



ih 



Ot 



Hx,t) 



2m 



{4 



+ (at; 



af): 



1 



-V 2 



$(x,t) 
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} 



(90) 



The case (af , ) = (1,1) (or equivalently (a^ , i 
and the case (af,a 



[1, 1)) recovers the usual linear Schroedinger equation, 



,. , - (0,1) recovers Eq. (1801) . 
Analogously, the nonlinear Klein-Gordon equation (j83j) can be further generalized as follows: 



V 2 $(x,t) 



1 d 2 ${x,t) 



+ ^^,t){af G + (af G -af G )g 



$(£, t) 



2(9-1) 



}■ 



(91) 



The case {af G ,a^ G ) = (1,1) (or equivalently (af G ,q) = (1,1)) recovers the usual linear Klein-Gordon 
equation, and the case (af G , af G ) = (0, 1) recovers Eq. 



Finally, the nonlinear Dirac equation 
d$(x,t 



ih 



can be further generalized as follows: 

D 



dt 



+ ihc{a-V)<S>(x,t) = /3mc 2 (%{af + (af - af )[^> j (x,t)/a j ] q ' 1 ^ ®(x,t). (92) 



The case (af , a q 



lJ " L> 1 — 1,1) (or equivalently (af , i 



(1,1)) recovers the usual linear Klein- Gordon equation, 



l ) a <j 



(0, 1) recovers Eq. 



and the case (a 

At the light of the situations that have been considered in the present Subsection, we may ask 
the following two points: What is the precise class of functions f{x) for which Eq. |75[ ) applies? (See 



Jauregui <fe Tsallisl . l2010al |) Is it possible, following along the lines of Eqs \72$ and \73\j or any other path, 
to find exact solutions of Eqs. (G2J), (OP and l9^i) ? 
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2.5. Dependence of q on the interaction-range of many-body Hamiltonians 

The BG entropy and its associated statistical mechanics are known to be very useful for many-body Hamil- 
tonians with say two-body interactions that are short-ranged and that do not introduce severe frustration 
(which would break down ergodicity, like it happens for instance for spin-glasses). A paradigmatic classical 
Hamiltonian system which violates (fo r < a < d) the sho r t-rang e-interacting condition (and introduces 
no frustration at all) is the a-XY one [Anteneodo Tsallisl . 1998^ | Fl: 



lo 9 1 1 — cos(6i — 9j) . „. . „, 

i=l i,j hJ 

where the N rotators are located at the sites of a d-dimensional simple hypercubic lattice (with periodic 
boundary conditions). The distances rij take the values 1, 2, 3, ... for d = 1, the values 1, v2, 2, ... for d = 2, 
the values 1, \/2, y/3, 2, ... for d = 3, and so on; due to the periodic boundary conditions, more than one 
value of rij can be defined between two given sites i and j (e.g., for d = 1, two such values exist generically) : 
in all cases we take into consideration only the smallest value of rij for a given couple N is defined 

as follows: 

j 



The role played by N is analyzed in [Anteneodo &: Tsallisl . [1998] . It is introduced here to conform to the 



vast literature existing for this Hamiltonian. It makes the Hamiltonian 71 n to be extensive for all admissible 
values of a. 

The (infinitely degenerate) fundamental state of the Hamiltonian (|93p corresponds to all rotators 
being parallel, and the corresponding total energy Un vanishes. We may define the asymptotic energy per 
particle u = limjv_ >00 Un/N. A critical value u c exists such that the system is ferromagnetically ordered 



for < u < u c , and it is 



paramagnetica lly disordered for u > u c . For example, for d = 1, it is u c = 3/4, Va 



(see [Campa et al l200ll : iTsallisl . T2009ai ] for further details). 



00 



For all values a/d > 1, only one collective thermal equilibrium exists, correctly (and analytically) de- 
scribed within BG statistical mechanics. For those systems, the limits limjv_ 5>00 lim^oo and lim^oo limjy 
of all thermostatistical properties commute. For all values < a/d < 1, two collective thermal station- 
ary states exist, namely the thermal equilibrium (corresponding to the ordering first t — > 00 and then 
— > 00), and the so-called quasi- stationary state (QSS) (corresponding to the ordering first N — > 00 and 
then t —> 00). Thermal equilibrium is correctly described within the BG theory. Not so the QSS, which 
emerge at values of u slightly below u c (e.g., for u ~ 0.69 for d = 1), and are not described correctly within 
the BG theory. Indeed, for QSS, the microscopic dynamics is nonergodic since the ensemble averages and 
the time averages do not coincide. The ensemble averages are neither BG (q = 1) nor q ^ 1 well described 
(approaches based on the Vlasov equation might be relevant). But the time a verages (by far the most 
relevant in experiments) appear to follow (/-statistics [Pluchino et a i l2007Ll200i |. Indeed, the distribution 



of momenta {pi} is definitively non Gaussian, and is well fitted by g„ e j-Gaussians with q w? > 1, where vel 
stand s for velocities. For t he a = mode l (known in the literature as the HMF model Antoni &: Ruffol . 



we have q ve i ~ 1.5 [Pluchino et al I2007I . l2008l |. In principle one expects q ve i(a,d) with q ve i(0,d) 



being independent from d. It is however, quite probable that q ve \ only depends on the ratio a/d, i.e., 
Qvei{ot/d) with g„ e ;(0) — 1-5. Moreover, it seems plausible that q ve i(a/d) monotonically decreases from 
Qvei(0) to q ve i(l) = 1 when a/d increases from zero to unity. 

A directly related question is the following. It is clear that, at the QSS, the corresponding (canonical) 
distribution for the entire phase space can not be given by the BG weight, i.e., we have that 

e -j8Wjv(pi,P2,.» >pjvA,02,~ ,0n) 

p( Pl ,p 2 ,..., PN ,0 1 ,9 2 ,...,e N ) + jj-j- JpjvdMe2j ^ r ,H ff(M2 • ( 95 ) 



9 Many other types of classical and quantum Ha miltonian systems can be thought of in l ong-range-interacting versions as 
illustrations along similar lines (see, for instance, [Nobre fc Tsallisl . 1 19951 : ICaride et all . ll983l p. 
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Indeed, if the BG weight was the correct one, we would have (since the kinetic and potential ener- 
gies commute) that the one- momentum marginal distribution P(pi) would be a Gaussian, namely the 



Maxwellian — 6 ■ But we numerically know that this is not true. Therefore the distribution 

f dpi e P p i< > 

P(pi,P2, — ,PNi #lj #2> —j &n) must be a different one. A possible Ansatz would be that, for N >> 1, 
it is given by 

-^stationary state n N(pi,P2,-,PN,ei,e 2 ,...,e N ) 

t~% / r\ r\ r\ \ *? stationary state /r\r>\ 

p(pi,P2,-,PN,0i,9 2 ,...,e N ) = ^ 57^7^ fl . a./. * ( 96 ) 



r , , , ,,, , n , a ~ ^stationary state U N (Pl,P2,- ■ ■ ,PN filfi 2, 

J dp 1 dp2...dp N d6 1 d6 2 ...d6 N e qstaUonary J ate 



where 



and 



{stationary state 



f(a/d) , (97) 



Stationary state = 9^/d, 0) . (98) 

Naturally we expect f(a/d) = 1 and g(a/d,/3) = f3 for a/d > 1. 

At this point, relevant open questions that arise are: Is it true that, at the QSS, the distribution in the 
full T phase space is given by Eq. ^96\) ? If so, what is the function q s tationary state = f(ce/d) for < a/d < 1? 
Is it true that the one-moment marginal distribution is a q ve i- Gaussian? Is it q ve i = q s tationary state, or, if 
not, what is it their relationship for < a/d < 1? Are q ve i(a/d) and q s tationary state{ a I d) universal in the 
sense that they are shared by a wide class of classical Hamiltonian systems? 

2.6. Temperature and the zeroth principle of thermodynamics 

The zeroth principle of thermodynamics plays a basic role in its axiomatic formulation. It states that if 
a system A is in thermal equilibrium with a system B, and B is in thermal equilibrium with a system C, 
then A is in thermal equilibrium with C. In other words, the concept of thermal equilibrium is transitive, 
which exhibits the great importance of the temperature, a quantity whose value is shared by all systems 
in thermal equilibrium. Can stationary (or quasi-stationary) states different from thermal equlibrium also 
satisfy this transitivity? It might be that some of them can. Indeed, in what concerns the behavior of two 
systems A and B at in QSS at somewhat different temperatures and being put in thermal contact suggests 
that. After contact between two equally sized systems (with N rotators and N(N — l)/2 links each), they 
evolve into a single double-sized system (with 2 N rotators) wh ose temperature is between the two initial 
temperatures: see Figs. H] and [5l and details in Tsallisl . 2009al |. All this occurs before the entire system 



makes the crossover to the BG regime, where thermal equilibrium takes place. 

The point which remains to be checked (for this and other similar long-ranged Hamiltonians) is: // we 
follow a full sequence of connections and disconnections between systems A, B and C, all of them at their 
respective QSS's, will the behavior be totally analogous to what is known to happen at thermal equilibrium 
(i.e., at the BG states)? 

2.7. Connection to thermodynamics and q-expectation values 

To generalize BG statistical mechanics for the canonical ensemble (from Tsallisl . 2009bl |). we optimize S q 



with the constraint 

w 

£> = 1 (") 

i=l 

and also 

w 

Y J PiE i = U q , (100) 
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System A 



System B 



Fig. 4. Systems A and B that will be put in thermal contact at a cer tain moment by allowing the coupling constant I to 
become different from zero. Here N = 5. From [Albuquerque et "al . |2008| . 




Time 



Fig. 5. Time evolution of the temperatures of A and B. The initial conditions are water bag for both A and B, at slightly 
different initial internal energies, hence slightly different initial temperatures. Her e N = 10000, and I is tak en zero until the 
moment indicated with a green vertical line, and I — 0.1 after that moment. From [Albuquerque et "al . 120081 ]. 



where 



P 



P 



IV 



i=i 



is the so-called escort distribution (Beck &: Schlogl . 19931 ]. It follows that p. 



p 1/q 



(101) 



There are various 



converging reasons for being appropriate to impose the energy constraint with the {Pi} instead of with the 
original {pi}. The full discussion of this delicate point is beyond the present scope. However, some of these 
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intertwined reasons are explored in [Tsallid . 12003 ]. By impo s ing Eq. (11001). we follow ITsallis et all . Il998 |. 
which in turn reformulates the results presented in |Tsalligl . Il988l : ICurado fc Tsallisl . \l99^. The passage 
from one to the other of the various existing for mulations of the above optimization problem are discussed 
in detail in [Tsallis et at Il998l : iFerri et aiEoO^ . 

The entropy optimization yields, for the stationary state, 



-fi q {Ei-U q ) 



Pi 



(102) 



with 



and 



A* 



ft 



T w v 



w 



(103) 



(104) 



(5 being the Lagrange parameter associated with the constraint (jlOOp , Eq. ()102j) makes explicit that the 
probability distribution is, for fixed /3 q , invariant with regard to the arbitrary choice of the zero of energies. 
The stationary state (or (meta)equilibrium) distribution (j!02|) can be rewritten as follows: 



with 



and 



Pi 



Z ' q 



Z ' q 



w 
i=i 



Pa 



Pa 



(105) 



(106) 



(107) 



' 1 + (1 - q)P q U q ' 

The form (| 105j) is particularly convenient for many applications where comparison with experimental or 



V(9-i) 



computational data is involved. Also, it makes clear that pi asymptotically decays like 1/E, 
and has a cutoff for q < 1, instead of the exponential decay with E\ for q = 1. 

The connection to thermodynamics is established in what follows. It can be proved that 

1 _ dS q 
T dU q ' 

with T = l/(k/3). Also we prove, for the free energy, 



for q > 1, 



(108) 



U„ ~ TS„ 



where 



ln g Z g 



\n q Z q 



■p^ q z q , 



(109) 



(110) 



This relation takes into account the trivial fact that, in contrast with what is usually done in BG statistics, 
the energies {Ei} are here referred to U q in (|1Q2[) . It can also be proved 

d 



as well as relations such as 



C a 



T dS q 



dU„ 



dT dT 



dT 2 



(111) 



(112) 
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In fact, the entire Legendre transformation structure of thermodynamics is g-invariant, which is both 
remarkable and welcome. 

Important questions that remain to be clarified include: What are the connections of quantities such as 
(3, f3 q , f3', and U q , with the thermo statistical quantities naturally appearing in models such as those focused 
on in Subsection \2. 51 ? 



2.8. q-triplet and possibly associated algebras 

For the most basic quantities (e.g., sensitivity to the initial conditions, relaxation towards equilibrium of 
correlation functions, equilibrium distribution of energies) of systems described by the BG theory, the ex- 
ponential function emerges ubiquitously as the adequate one. This function is replaced by the ^-exponential 
one for systems described by nonextensive statistical mechanics. The question appears about what is the 
value of q to be used. A wide number of examples show that the value of q is directly associated with the class 
of properties that are being studied. For example, for dissipative one-dimensional maps at the edge of chaos 
we have that the sensitivity to the initial conditions is characterized by qsensitivity 

< 1, the entropy produc- 
tion is characterized by q en tropy production (and for such systems we verify that q se nsitivit y = Qentropy production), 
the relaxation of the entropy towards its saturation value is characterized by q re iaxation > 1> the sums of 
many successive iterations are characterized by q a ttractor > 1, and so on. In other words, for a given 
nonextensive system, we typically have not one value of q, but an infinite number of them, most probably 
interconnected in such a way that at the end only one (or very few) are independent, thus characterizing 
universality classes of nonextensivity. In the limit when a BG regime is approached, all these values of q 
are typically expected to merge into the single value q = 1. 

Very little is known nowadays about the (plausible) analytical con nections of such indices with the 
following algebra emerging within the g-CLT [Umarov et all . 20081 . 2010 ] and related matters: 



! I m 

+ — (0 < a < 2; m = 0,±1,±2,±3, ...) , (113) 



1 — q i 1 - 9o « 



and 



+ — (0 < a < 2; m = 0,±1,±2,±3, ...) , (114) 



1 — Q i Qo — 1 a 

where a = 2 corresponds to g-Gaussians, and < a < 2 corresponds to (q, a)-stable distributions. See 
Fig. [6J As examples of g-triplets and analogous relations, let us mention, among others available in the 
literature, the following ones: 



2005]): 



i) Fluctuations of the magnetic field of the solar wind (as detected by Voyager 1 [Burlaga Vinas 



Qsensitivity — 0.6 it 0.2 , 
Qreiaxation — 3.8 i 0.3 , 
Qstationary state = 1.75 ±0.06, (115) 

possibly interpreted Tsallis et all . 20051 ] as 

Qsensitivity — 1/2 , 
Qreiaxation 4 , 
Qstationary state = 

7/4. (116) 
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2 3 4 5 



V.v, 




Fig. 6. Curves of <?^/ Q (top) and 1 m / a (bottom) as functions of go for typical values of m/a. They are respectively given 
by Eqs. (|113[) and (|114[) . The three black dots in the top figure correspond, from top to bottom, to q^ = q re laxation = 4, 
'stationary state = 7/4 and q£ = q sensit ivity = -1/2 respectively [Tsallis et a 1, l2005| ]. 



4 
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(ii) Edge of chaos (Feigenbaum point) of the logistic map (see Tirnakli et all . l2010l ; lFuentes &; Robledol . 



2010] and references therein): 



Qsensitivity 
^relaxation 
Istationary state 



0.244487701341282066198... 

2.249784109..., 

1.65 ±0.05. 



(117) 



(iii) Fluctuations of the width (above Buenos Aires) of the Ozone layer |Ferri et"al . 2010( | : 



Qsensitivity — 8.1 zb 0.2 , 
Qrelaxation = 1-32 i 0.06 , 
Qstationary state = 1-89 ± 0.02 . (H8) 



As we can see, these examples suggest that quite widely qsensitivity < 1 < qstationary state < qrelaxation- 
As open questions we may emphasize: What are the precise connections between physical properties 
and the q-triplet and similar quantities? How many of those indices are independent? Are there direct 
connections between these indices and the q-CLT algebra shown in Eqs. M1'J\) and \11J$ ? 



2.9. Universality classes, classification of entropies 

The BG entropy and its associated exponential distribution for thermal equilibrium have been extended, 
during the last two decades, in the sense of thermodynamics and statistical me chanics , into other forms. 
Stati onary states have been discussed which cor rep ond to the a - exponenti al form iTsallis . 1988ll. logarith mic 
form [Curadol . 1999 : Curado fe Nobrel. 12004 1 , stretched- exp o nential for m Anteneodo &: Plastinol . Tl999| . as 
well as other, more general, forms | H aneifc Thurnerl. Eoill: ff^Q . M~ 

The present belief is that BG entropy and statistics are sufficient but not necessary for thermodynamics. 
In other words, thermodynamics might be more powerful than the role attributed to it by BG statistical 
mechanics. A question which arises then naturally is: What is the most general form of entropy which can 
be consistent with thermodyna mics, more precisely with the zeroth, first, sec ond, and third principles? What 
are the super statistical forms Beck h Cohen , 20031 : Tsallis fc Souza , 2003] which would correspond to this 
general entropy and statistics? 



3. Conclusions 

During the last two decades the nonadditive entropy S q and its associated nonextensive statistical mechanics 
have been an d are being intensiv ely studied world wide (over three thousands of papers are available in 
the literature Bibliography . 2010l |). A variety of analytical, computational, experimental and observational 
results and applications provide nowadays what one may consider as a relatively clear understanding of 
its domain of validity. Nevertheless, various relevant points still remain to be clarified. We have presented 
here a set of such open questions in the hope that future work will improve our insights. 
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